Direct bandgap optical transitions in Si nanocrystals 
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The effect of quantum confinement on the direct bandgap of spherical Si nanocrystals has been 
modelled theoretically. We conclude that the energy of the direct bandgap at the T-point de- 
creases with size reduction: quantum confinement enhances radiative recombination across the direct 
bandgap and introduces its "red" shift for smaller grains. We postulate to identify the frequently 
reported efficient blue emission (F-band) from Si nanocrystals with this zero-phonon recombina- 
tion. In a dedicated experiment, we confirm the "red" shift of the F-band, supporting the proposed 
identification. 

PACS numbers: 73.22.Dj, 78.67.Bf, 78.67.Hc 
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Crystalline silicon, which dominates electronic and 
photovoltaic industry, has notoriously inferior optical 
properties for light emission. This follows from its en- 
ergy structure which features an indirect bandgap, with 
the absolute minimum of the conduction band being very 
much displaced — Ak = 0.85 x [2% /a) along the [100]- 
direction — from the T-point, where the maximum of 
the valence band is positioned. The direct bandgap en- 
ergy at the T-point amounts to E^" — 3.32 eV in bulk 
Si. On the way towards photonic and opto-electronic ap- 
plications of silicon, possibly the most promising route 
is offered by space-confinement-induced changes of en- 
ergy structure. Following the pioneering paper of Can- 
ham on porous Si [l[, investigations of various forms of 
nanostructured Si have been undertaken. In particular, 
very interesting results were obtained for Si nanocrystals 
(NCs) where strong effects of quantum confinement have 
been observed for grain sizes comparable to or smaller 
than the effective Bohr radius of 4.3 nm — e.g. see 
Ref. for an extensive review. It has been shown that 
while the energy structure in Si NCs retains its indirect 
character, quantum confinement leads to relaxation of 
the momentum conservation requirement, thus making 
zero-phonon optical interband transitions possible. Nev- 
ertheless, phonon-assisted radiative transitions dominate 
for nanocrystal sizes down to 2 nm ||. Therefore, the 
emission from Si-NCs has two characteristic features: for 
smaller grains, the spectrum shifts to the blue and its in- 
tensity increases. In the literature, this widely tuneable 
and rather efficient emission due to ground state electron- 
hole recombination is usually referred to as the "slow" 
band (5-band). It has attracted a lot of attention and 
prospects of Si-NCs- based laser have been discussed 0] • 

In addition to the 5-band whose origin is well estab- 
lished, optical investigations of Si-NCs revealed also an- 
other emission band at higher energies. This "blue" 



band is characterized by much faster decay dynamics 
(hence the name -F-band), ranging from pico- to nanosec- 
onds. A similar emission band has also been reported for 
porous Si [5j. Its origin is currently debated both ex- 
perimentally and theoretically ■ Usually, it is postu- 
lated to arise due to radiative recombination of electrons 
localized at oxygen-related states at Si/SiC>2 interface 
[j| [l(| • This interpretation has been recently challenged 
by Valenta [ll[ , who presented evidence that the F-band 
should rather be identified with transitions related to the 
NCs themselves. Emission bands in the visible have also 
been reported in colloidal Si-NCs [HI, [l3[ and postulated 
to arise due to direct transitions in view of their short 
lifetime. 

In this paper, we consider effect of quantum confine- 
ment on the conduction band at the T-point in spherically 
symmetric Si-NCs. According to the presented theoreti- 
cal model, quantum confinement applied to the conduc- 
tion band states in the vicinity of the T-point, creates two 
series of states with energies higher and lower than the 
energy of the conduction band at the T-point in bulk Si 
(3.32 eV), forming a gap between them. We attribute the 
lower series of states as the origin of direct optical transi- 
tions responsible for the F-band. In that interpretation, 
the F-band is microscopically identified as arising due to 
zero-phonon recombination of "hot" electrons from the 
T-point of the conduction band. We discuss that as the 
confinement gets stronger with decreasing the NC size 
and the direct band gap shrinks, the "red" shift of the 
F-band appears. We confirm this theoretical conclusion 
by dedicated experimental evidence concerning spectral 
dependence and dynamics of the F-band in a dense solid 
state dispersion of Si-NCs in a Si02 matrix. We demon- 
strate that the proposed microscopic identification of the 
F-band is consistent with this experimental data as well 
as with the available literature. 
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FIG. 1: Fig. 1. Right panel: schematical representation of the 
dispersion relations for the conduction and the valence bands 
of bulk silicon along the [001] direction [TJ]- Left panel: quan- 
tum potentials in real space and confinement energy levels for 
electrons and holes. Our model attributes the lower quantum 
levels (solid lines) of electrons to the states originating from 
A-valley of the conduction band, and the higher excited ones 
(dotted lines) to the states of the F-point. 



The scheme of the lowest conduction band in bulk sil- 
icon (based on detailed calculations from Ref. [14|) is 
shown in the right panel of Fig. 1. In the center of 
Brillouin zone its states are of Tis symmetry, thus be- 
ing 3- fold degenerate (without spin taken into account), 
and they are related to two subbands, one of which 
(the "heavy" one) has negative effective mass. It is the 
negative effective mass that leads to emerging of space- 
confinement levels below the band edge and therefore to 
decrease of the threshold of direct optical transitions. 

Counting energy of electrons from the Tis-point, one 
can write a generalization of the Luttinger Hamiltonian 
in spherical approximation: 



where too is the free electron mass, p is the momen- 
tum operator, and J is the unitary angular momentum 
operator acting in the space of Bloch amplitudes. The 
coefficients A and B, related to the effective mass values, 
are extracted from the results of Ref. [l4j by weighten- 
ing the band parameters obtained for different directions: 
A = 0.4^ooi + O.6A111 = 0.58, B = 0AB ool + 0.6B m = 
—0.85. Such values correspond to positive energy and 
"light" effective mass of m e i — 0.35 to for doubly degen- 
erate subband and negative energy and "heavy" effective 
mass of m e h = —0.45 mo for the other nondegenerate 
subband. In a spherical quantum dot Hamiltonian ([1]) 
results in three types of states [l6| , which are the eigen- 
functions of the square F 2 of the full angular momentum 
operator F = J + L (operator L = — ir x V acts on the 




NC diameter (nm) 

FIG. 2: Fig. 2. Energy levels of electrons at the T-point 
in spherical Si-NC as a function of its diameter. According 
to Fig. 1, zero energy corresponds to the energy of bulk Si 
conduction band at the T-point. 

envelope parts of the wave functions only) and its pro- 
jection F z . They are 

i) "heavy" states Neho with F = 0, 

ii) "light" states Nel F (F > 1), 

iii) "mixed" states ±Nem F (F > 1), 

where N is the main quantum number and the subscript 
shows the value of F. 

To investigate qualitatively the effect of confine- 
ment on the energy spectrum, we calculated the space- 
quantization levels of electrons in the vicinity of the Pin- 
point assuming for them infinitely high energy barriers 
at the NC boundary. The resulting energy spectrum as 
a function of NC size is shown in Fig. 2 . One can see 
that the spectrum is split: "light" states have positive en- 
ergy, "heavy" — negative, and "mixed" ones have both 
positive and negative energy levels. The most important 
feature of the resulting spectrum is that the energy gap is 
formed and the absolute values of negative energy levels 
are rising with decreasing diameter of NCs. We should 
remark that the appearance of such a gap was predicted 
in Ref. [15| with help of rudimentary pseudopotential cal- 
culations. The quantization energies of holes are rising 
as well [lH, but it does not change the trend, and the 
transition energy does decrease for given transition with 
decreasing the size. 

The rate of direct spontaneous optical recombination 
is given by [171 ]: 

r'd = ^a 3 Lm ont \M\ 2 F 2 , (2) 

where a = h/(moaoc) rs 1/137 is the fine structure con- 
stant, ao is Bohr radius, c is the speed of light, fku is 
the transition e nerg y, T = 3nl ut /(nf n + 2n 2 nt ) is the lo- 
cal field factor [17j, with n- ln and n out being refraction 
indexes of media inside and outside the NC respectively 
(for the NCs under consideration T ~ 0.42). The dimen- 
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TABLE I: Factors l\ for several interband transitions, as 
well as their energies for two NC sizes: 3 and 7 mil. 







d = 


3 nm 


d — 7 nm 


Transition 


T 2 


hw 


1% hw 


le/io — 


> 1/W712 


1.185 


2.95 eV 


1.277 3.28 eV 


— \em\ - 


->• lft/l2 


0.099 


2.53 eV 


0.130 3.23 eV 


— lemi - 


■¥ lhmi 


0.016 


2.17 eV 


0.002 3.13 eV 


2eh - 


> \hu12 


0.077 


1.45 eV 


0.047 3.01 eV 



sionless matrix element M in ([2]) is defined by: 



IMF 



h 2 N c . N h 



E 

M B ,M h 



E 

a=x,y,z 



E / (^T^d S r{u v \p a \u c ) 



> (3) 



where N = 2F C + 1 and N h = 2F h + 1 arc the 
degeneracy factors of the electron and hole states re- 
spectively, and 7/>* fo ( -i /'^ fh ) are the envelope parts of 
the electron (hole) wave function \& f c m c = J2 c " l i J c I " u ' c 
(^F h M h = ^ytp^Uv)- The Bloch amplitudes basis u c 
is chosen in the form of spherical components uq = Z c , 
u± = =F^/l/2(X c ± iic), where AT C = x, Y c = y, and 
Z c = z are the functions of the T 15 representation. It 
should be noted that the consideration of hole states 
in [l6j has also been produced with the Bloch ampli- 
tude basis u v in the form of spherical components but 
using functions X v — yz, Y v = xz, and Z v — xy of the 
T25' representation. The difference in the Bloch ampli- 
tude bases makes the direct optical transitions between 
electron and hole states permitted. 

The value of the matrix element of momentum between 
the Bloch amplitudes of the valence and the conduction 
bands p cv allowed by the selection rules can be writ- 
ten in the form: p 2 cv = \{X c \p y \Z v )\ 2 = \(Y c \p z \X v }\ 2 = 
\{Z c \p x \Y v )\ 2 = Q 2 h 2 /al where Q 
total matrix element ([3]) is \M\ 2 



1.07 [14|. And the 

Q 2 Ir, where Ir is 



a combination of overlap integrals of the envelope func- 
tions. There are also selection rules for overlap integrals 
of envelope functions in Eq. ([3]) arising due to the spher- 
ical symmetry of both electron and hole states. E.g., 
in our consideration, radiative transitions from "heavy" 
electron states Neho can only go into the valence band 
states of "mixed" type, characterized by F — 2: Nhni2- 
These rules should not be strict for real NCs, which are 
not perfectly spherical, and transitions disabled by our 
model are still possible, but we expect them to be much 
weaker. 

Finally, for the estimation of radiative recombination 
rate one can use: 



= 1.0 x 10 8 



3cV 



Ir, s- 



(4) 



The values of I R for the radiative transitions with the 
smallest energies are listed in Table Q] For the most ef- 



fective transition Neho Nhm,2 the radiative life time is 
of the order of 10 ns. We propose that the emissions ob- 
served as the F-band arise due to transitions in Table Q] 
Table Q] also illustrates that the energy range of possible 
transitions increases with decrease of NC diameter. 

The radiative emission under consideration can be pro- 
duced only by "hot" confined electrons. Thus the rate of 
nonradiative relaxation of these electrons is a key issue 
for the observation of this emission. We have studied 
transitions between levels in the vicinity of the Tis-point 
taking place due to multiple emission of optical phonons 
and emission/absorption of one promoting acoustic lon- 
gitudinal phonon. The probability of such a transition is 
given by W — r^ 1 J p (T, S), where r ac is the time deter- 
mined by the interaction with an acoustic phonon, T is 
temperature, S is the Huang-Rhys factor, and J P (T,S) 
is the factor specific for multiphonon transitions which 
for S <C 1 and kT less than the optical phonon en- 
ergy huj ~ 60 meV can be well approximated by S p hoi , 
where p is the number of emitted optical phonons |18| . 
We have calculated the Huang-Rhys factor following 
Ref. [H[ and using Bir-Pikus Hamiltonian i? e °Pp h (u opt ) = 

71 V K pt {Jrh} + u$%J z J,} + for the 

interaction with optical phonons, where 2{J Q Jg} = 
J a J/3 + JpJa, a is the lattice constant, rf = —16.9 eV [l9| 
is the interaction constant, and u opt is the relative atomic 
displacement induced by the optical phonon mode. In re- 
sult we have S = 0.13/c? 3 for le/io — > — lemi transition, 
S = 0.04/rf 3 for -lemi -> -lem 2 , and S = 0.22/d 3 
for — lem2 ^ 2eho, where the NC diameter d should 
be taken in nm. Further we have calculated t^} and 
then W taking into account the energy difference which 
is not compensated by the optical phonons and using 
the interaction Hamiltonian H^ ph = a c Vu ac , where 
a c = —10 eV [lj| and u ac is the atomic displacement 
induced by the acoustic phonon mode. The highest tran- 
sition rates are of the order of 10 9 s _1 for the largest NCs 
considered here. They decrease rapidly with the NC size. 
One can notice that for transitions letiQ — > — lemi and 

— lemi — > — lem.2 several optical phonons are needed even 
for large NCs with extremely small Huang-Rhys factors 
(see Fig. 2). In case of small NCs more phonons are 
required and Huang-Rhys factors are still small so that 
the phonon-induced relaxation is suppressed. The state 

— lemi can decay quicker than the other states due to 
smaller energy spacing to the lower neighboring state but 
it anyway does not contribute to the considered radiative 
transitions (see Table I). 

The proposed theoretical identification of the F-band 
is directly supported by photoluminescence (PL) experi- 
ments performed on Si-NCs embedded in a Si02 matrix. 
Two different samples were used, both prepared by ra- 
dio frequency co-sputtering and subsequent annealing, 
resulting in an average NC diameter of 3 nm (sample 
A) and 5.5 nm (sample B), with a small size distribu- 
tion — see, e.g. Refs. [10, [2l[ for more details on sample 
preparation and measuring techniques. Figure 3 shows 
the results obtained for both samples in two differently 
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FIG. 3: Fig. 3. (a) S (scattered symbols) and F(solid line 
+ symbols) PL bands for sample A (circles) and sample 
B (squares). Dashed arrow indicates the position of direct 
bandgap recombination as experimentally measured for bulk 
Si [2J]. (b,c) PL dynamics for sample A (figure b) and sam- 
ple B (figure c): S-band decay is displayed in main pan- 
els. The kinetics were fitted with a stretched exponential: 
IpL=Ioe" (t/T) ' 3 , with = 0.8 The relevant time con- 

stants are rf « 20 (is for sample A and r| « 5f0 fxs for 
sample B. The inset shows dynamics of the F-band in a time 
window of 50 ns taken at the maximum intensity values of 
2.5 eV and 3 eV for samples A and B, respectively. 

conducted PL experiments. With the steady-state tech- 
nique, the slow (microsecond range) time-integrated PL 
signal was monitored. This is dominated by the low- 
energy 5-band (indicated in Fig. 3(a) with O ) and, as 
a result of quantum confinement, shows the well-known 



shift towards higher photon energies for smaller NCs. 
The spectral profile of the F-band has been obtained by 
taking the initial amplitude of the ps-resolved dynamics 
and is illustrated in Fig. 3(a) ( — ). As can be seen, the 
F -band shows a red shift upon size decrease — an effect 
opposite to the blue shift characteristic for the 5-band. 
This behavior is in agreement with the presently calcu- 
lated shift towards lower energies of the electron energy 
levels at the T-point — see Fig. 2 — and the consequent 
reduction of the direct bandgap value. We note that 
the linewidth of the F-band clearly increases for smaller 
NCs, as indeed expected due to quantum confinement. 
The time-resolved PL measurements have been obtained 
in a time-correlated single photon counting experiment. 
Figs 3(b) and 3(c) illustrate dynamics of S and F PL 
bands for both samples. The decay time of the 5-band 
(main panels, nanosecond resolution) is of an order of 
~ 10 2 fj,s, i.e. similar to the estimated radiative re- 
combination time of excitons in Si-NCs 0, [22J. It is 
somewhat faster for smaller NCs (sample A in respect 
to B], reflecting the increased radiative recombination 
rate [23|]. The decay dynamics of the F-band (insets, pi- 
cosecond resolution) is characterized by a time constant 
of rp ~ 500 ps, similar for both samples. Finally, we 
point out that the fact that the F-band can be monitored 
using (ps-resolution) time-resolved PL experiments, im- 
plies that an upper limit can be set to the ratio of the 
non-radiative transition rate to the radiative recombina- 
tion rate of (r^Jy 1 (r^)" 1 < 10 2 . On that basis, 
the radiative recombination time for the F-band can be 
estimated as r r ^ d < 50 ns. 

We conclude that the characteristic features following 
from the proposed theoretical model for the direct T- 
point phonon-less recombinations — (i) decrease of the 
energy and (ii) linewidth broadening for smaller NC sizes 
are experimentally confirmed for the F-band emission 
from Si-NCs. Also the radiative lifetime of the F-band 
estimated from the measured PL dynamics agrees with 
the calculated values. These additional new experimen- 
tal evidence taken together with the information on F- 
band available in the literature, agree very well with the 
proposed identification of this fast emission with radia- 
tive recombination across the direct bandgap of Si, whose 
energy has been lowered in Si-NCs due to quantum con- 
finement. 

Further research, currently on the way, will elucidate 
details of the direct bandgap tuning in Si-NCs. 
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